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Abstract 

Two-dimensional scattering by homogeneous and layered dielectric elliptical cylinders is analyzed 
following an analytical approach using Mathieu functions. Closed- form relations for the expansion 
coefficients of the resulting electric field in the vicinity of the scatterer are provided. Numerical 
examples show the focalizing effect of dielectric elliptical cylinders illuminated normally to the axis. 
The influence of the confocal dielectric cover on the resulting scattered field is envisaged. 

PACS numbers: 41.20.Jb, 42.25.Fx, 42.25. Gy 



* Electronic address: ecojocaru@theory.nipne.ro 

I 







I. INTRODUCTION 



The problem of wave scattering by an elliptical cylinder has been known for some time 
[U [2] . Exact solutions can be obtained by separating the wave equation in elliptic cylindrical 
coordinates and constructing exact eigenfunction expansions. But these exact solutions have 
been of limited practical value due to problems associated with the computation of the 
corresponding eigenfunctions, i.e. angular and radial Mathieu functions. 

Recently, a new efficient method for computing these functions was presented [3j, the 
method being implemented in matlab j4j. To extend the applicability of the new computa- 
tional method, we derive in this paper the eigenfunctions for two-dimensional (2D) scattering 
by dielectric elliptical cylinders. Closed-form relations for the normally incident plane wave 
scattering by homogeneous and layered elliptical cylinders are provided. 



II. WAVE EQUATION IN ELLIPTIC CYLINDRICAL COORDINATES 

Firstly, the elliptic cylindrical coordinates (u,v,z), as shown in Fig. [TJ are introduced in 
terms of rectangular coordinates (x, y, z) 

x = /cosh u cos y = fs'mhu smv, z = z, (1) 

with < u < oo, < v < 27T, and / the semifocal length of the ellipse. The contours 
of constant u are confocal ellipses, and those of constant v are confocal hyperbolas. The 
z axis coincides with the cylinder axis. In the following we consider the transverse-electric 
(TE) polarized electromagnetic field (i.e., the electric field only exists in the z direction). 
By expressing the Laplacian in elliptic cylindrical coordinates pQ, the scalar wave equation 
in a homogeneous dielectric medium of permittivity e is given by 

2 fd 2 E z + d 2 E z \ , / 1 d 2 E z d 2 E z 



fiof 2 (cosh2u — cos2v) \ du 2 dv 2 ) Ve /x dz 2 dt 2 

An exp(-iujt) time dependence is assumed, where uj is the circular frequency. Separation of 
variables implies that we assume a solution E z of the form 

E z = Z(z)S(v)R(u) (3) 
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On substituting Eq. ^ in Eq. ^ we obtain 



Id? 

A 2 

dv 
dt£ 



l + fc*)z(z) = 0, (4a) 

2 + (a - 2g cos 2v) ] S(v) = 0, (4b) 
2 - (a - 2<? cosh 2u) i?(^) = 0, (4c) 



where k z and a are separation constants, q = f 2 k 2 /4, with k 2 = k 2 — k 2 z) where k = k y/e, 
k = 27r/A, and A is the wavelength of the incident light in vacuum. By taking the advantage 
of the elliptical cylinder geometry, we focus our analysis on the 2D scattering of normally 



incident plane waves (k z = 0). Equation (4b) is known as the angular Mathieu equation. 



The solution is denoted by S pm (q, v, n) where p, m denote even (e) or odd (o), and n denotes 



the order [3j 0j. Equation (4c) is known as the radial Mathieu equation. Similarly to the 
circular cylindrical coordinates where the radial solution is expressed in terms of Bessel 
functions J n (-), Y n (-), Hn\ m ), and Hn\-), the radial Mathieu equation has four kinds of 
solutions: J pm , Y pm , H pmU and H pm2 , where p, m = e, o [D El 0]. 

Let a scalar plane wave of unit amplitude be normally incident on the axis of the elliptical 
cylinder in air, so that the propagation direction forms an angle (j) with the x-axis in the 
x — y plane. Hence the electric field of the incident plane wave is given by 

jjjin ^iko (x cos 4>+y sin <f>) ^iko f (cosh u cos v cos 0+sinh u sin v sin cf>) 

To obtain exact solutions for the scattered wave, we expand the incident field in Eq. ^ in 
terms of the eigenfunctions of the elliptical cylinder [TJ [3] 

E™ = V8n ^2 inj P m(q, u, n)S pm (q, v, n)S pm (q, 0, n)/N pm (q, n), (6) 

n 

where N pm is normalization constant. 

III. SCATTERING BY HOMOGENEOUS ELLIPTICAL CYLINDERS 

Consider a homogeneous dielectric elliptical cylinder of permittivity e\ with the boundary 
located at u\, in air. The electric field E z outside and inside the elliptical cylinder is given 
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by [51 E] 

(u > Mi) E z = V8ny~] i n J pm (qo, u, n)S pm (q , v, n)S pm (q , (p, n)/N pm (q , n) 

n 

+ i n a p s £>(n)H pml (q , u,n)S pm (q ,v, n)S pm (q , (j), n), 
(u<ui) E z = V^^i n a p 1 l(n)J pm (q 1 ,u,n)S pm (q 1 ,v,n)S pm (q ,(j),n), (7) 

n 

where Opli and a p s ^ are coefficients to be determined from the boundary conditions. The 
tangential magnetic field H v in each region is obtained from E z [U EJ [3], 

(u>u 1 ) H v = j^[-tueo^ 

1 / dE \ 

(u< Ul ) H v = w - p (-^u J e 1 ^), (8) 

where p = /(sinh 2 2i + sin 2 ?;) 1 / 2 . The boundary conditions require the continuity of the 
tangential components of the electric and magnetic field at the boundary u — u\. One 
obtains 

yZ in i J prn(qo, u u n)/N pm (q , n) + a { p s ^(n)H pml (q , u u n)]S prn (q , v, n)S pm (q , (f), n) 

n 

= ^2 i n Upm{n)Jprn{qi,u u n)S pm (q u v, n)S pm (q , (j), n), 

n 

^2 in [ J pm(<?o, Mi, n)/N pm (q , n) + a p s £(n)H' pml (q , u u n)]S pm (q , v, n)S pm (q , <f>, n) 

n 

= ^ ' n «im( n )^m(?i. U U n)S pm (qi, v, n)S pm (q , (f), n), (9) 

n 

where the prime denotes differentiation with respect to u. To determine the coefficients 
4mH and a p S m (n), m = e, o, we multiply Eq. ^ by S pm '(qo, v, n), integrate the resulting 
equation over v from to 27r, and use the orthogonality relations of the angular Mathieu 
functions (HEJHj. Applying those orthogonality relations allows waves in each order n, and 
each combination p, m = e, o to decouple. One obtains 

•/pro (go, ui,ri)/N pm (q , n) + a { p ^(n)H prnl (q , u u n) 

= ^pm( n hpm(n)J p m(qi,U U Tl)/N prn (q , 77,), 

4m u u n)/N pm (q , n) + a^(n)H pml (q , u u n) 

= ai pl( n ) r ypm{n)J'(q 1 ,u ll n)/N pm (q , rc), (10) 
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where correlation factors 7 pm (n), p ) m = e, o are defined by 

/>2tt 

/ 5 pm /((?o,^,n)5' pm ((?i,^,n)d^ = 7 pm (rc)5 mm /. (11) 
Jo 

Expressions for the correlation factors in terms of expansion coefficients of Mathieu functions 



are given in [4]. From Eq. (10) we find 

(i) / x _ Jpm(qo,v>i,ri) J V m(%^ x ,n) - Jtfpmi(qo,Ui,ri) 
pm l P m(n)J pm (q u u u n) / pm {(l\,n ll n) - J^ pm i(qo,u u n)' 

pm N prn (qo,n)H pml (q ,u u ri) f P m{q\^n) - J^ pm i(qo,u u n)' 

where we introduced the log-derivative function & = F "/ F \ F being J pm or H pml: the 



prime denoting differentiation with respect to u. 

As an example we consider an elliptical cylinder of permittivity e = 2 [Figs. [2](a) and (b)] 
and e = 3 [Figs. [2](c) and (d)]. The semifocal length is / = 0.1m. The elliptical boundary 
is at U\ = 1.4436 (shown as black contour in the figures) with semiaxes d x = 0.2236m and 
d y = 0.1414m. The scatterer is illuminated from — x direction, perpendicularly to axis z. 
The incident light wavelength is A = 0.15m. We found that n = 7 is enough to assure 
the convergence of the calculated scattering field [4j. Snapshots of the normalized electric 
field [Realms)] distribution in the vicinity of the scatterer are shown in Figs. [2](a) and (c) 
for e = 2 and 3, respectively, normalization being against the maximum absolute real value. 
Figures [2](b) and (d) show the distribution of \E Z \ : normalized to the maximum value, for 
e = 2 and 3, respectively. One can see the focalizing effect of the perpendicularly illuminated 
dielectric elliptical cylinder, the greater the permittivity value, the higher the intensity on 
the focalizing direction. 



IV. SCATTERING BY LAYERED ELLIPTICAL CYLINDERS 

Consider a homogeneous dielectric elliptical cylinder of permittivity e\ and boundary at 
Ui covered by a confocal dielectric elliptical layer of permittivity e 2 and outer boundary 
at u 2) in air. Let q\ and q 2 be the respective elliptical parameters, q\ = pk^ei/A and 
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12 = / 2 ^o e 2/4- The electric field E z in each region is given by 
(u > u 2 ) E z = \/8% ^2i n [J P m(qo,u,n)/N pm (q ,n) 

n 

+ «pm( n )#pmi(<7o, u, n)]S pm (q , v, n)S pm (q , <j>, n), 
(ut<u< u 2 ) E z = V^^i n la p 2 l(n)J pm (q 2 ,u,n) (13) 

n 

+ ^pi(n)H pml (q 2 , u, n)]S pm (q 2 , v, n)S pm (q , (f>, n), 



(u < Mi) E z = V^^i n a p 1 l(n)J pm (q 1 ,u,n)S pm (q 1 ,v,n)S pm (q ,(f),n). 



1 / nW \ 

a =^( -*■*'-*:)■ }=1 - 2 - (14) 



The tangential magnetic field H v in each layer is 

_ luje . d ^l 
f>f,< il> V 3 du 

The boundary conditions require the continuity of the tangential components of the electric 
and magnetic field at the boundaries u = Uj, j = 1, 2. Following the same procedure as in the 
previous section we obtain four equations for each order n. (For simplicity, the dependence 
on order n is skipped.) 

J P m(qo,u 2 )/N pm (q ) + a p s ^H pml (q ,u 2 ) 

= 7pi/N pm (qo)[a p 2 lJ pm (q2,u 2 ) + a { p lH pm i(q 2 ,u 2 )}, 

Jp m {<lo,u 2 )/N pm (q ) + a p s ^H' pml (q Q ,u 2 ) (15) 
= 7pl/ N P m(qo)[a { p lJpm(Q2, u 2 ) + af m H' pml (q 2 , u 2 )}, 

l { piWp 2) m Jprn{q2,u 1 ) + a p 3 lH pml (q 2 , Ul )} = a$,7pm-Wgi,Mi), 

7S[a?i4m(92,«i) + a p 3 lH pml (q 2 , Ul )] = a^S^m^i, 

where correlation factors 7pm, p,m = e, o, j = 1, 2, are defined by 

/>2tt 

/ S pm >(q ,v)S pm (q j ,v)dv = 'y^S m m', J = 1,2. (16) 
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Relations (15) can be written in the form A • X = B, where 

/0 r rpmJ P m(Q2,u 2 ) Jp 2 iH pml (q 2 ,u 2 ) -H pml (q ,u 2 )N pm (q )\ 

lpmJ'{q2,U 2 ) 7pmflJ m i(</2,«2) ~ Km l(Q0, U 2 )N pm (q 



A = 







ui) 7p m H pml (q 2l Ui) 

\-JpmJ p m(qi, U l) IpmJ'pmi^Ux) 7pm#pml(«2, «l) 








/ 



X 



(2) 
pm 

(3) 
pm 

\ (sc) | 

\or P mJ 



a 



a 



and 5 



/ ^m(<?0,^2)\ 






(17) 



The unknown coefficients a4 J m (j = 1, 2, 3, sc) are determined at each order n by solving the 
matricial equation. In matlab, an equation of this form is solved simply with command 
X = A\B. 

As an example, we show the distribution of \E Z \ normalized to the maximum value for a 
layered elliptical cylinder with permittivities e\ = 2, e 2 = 3 in Fig. [3](a), and e\ = 3, e 2 = 2 
in Fig. [3](b). The semifocal length is / = 0.1m. The inner boundary is at u\ = 0.8814 
(with semiaxes d x = 0.2m, d y = 0.1m) and the outer boundary is at u 2 = 1.4436 (with 
semiaxes d x = 0.2236m, d y = 0.1414m), the boundaries being shown as black contours in 
the figure. The wavelength of the incident light is A = 0.15m, the light being incident 
from the — x direction, normally to the z axis. One can see that a confocal cover influences 
on the scattered field distribution; a lower (higher) permittivity cover decreases (increases) 
the intensity of the scattered field. The focalizing effect of the perpendicularly illuminated 
elliptical cylinder is increased by a lower permittivity confocal cover. 



V. SUMMARY 



In this contribution we have studied the plane wave scattering by homogeneous and 
layered elliptical dielectric cylinders following an analytical approach that is based on Math- 
ieu functions. By taking the advantage of the elliptical cylinder geometry, we focused our 
analysis on the 2D scattering of normally incident plane waves. We have provided closed- 
form relations for the expansion coefficients of the resulting electric field in the vicinity of 
the homogeneous and layered dielectric elliptical cylinders. Numerical examples have been 
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provided showing the focalizing effect of a homogeneous elliptical cylinder illuminated per- 
pendicularly to the axis, and the influence of a confocal cover layer on the resulting electric 
field distribution. The results presented could be a valuable contribution to the enlargement 
of the Mathieu functions applicability. 
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FIG. 1: The elliptic cylindrical coordinates. F\ and F2 are the foci of the ellipse; / is the semifocal 
length. 
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FIG. 2: (a) Snapshot and (b) norm of the resulting electric field distribution in the vicinity of an 
elliptical cylinder with permittivity e = 2; (c) and (d), the same with permittivity e = 3. The light 
is incident from the — x direction, in the x — y plane, normally to the z axis. 
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FIG. 3: (a) Norm of the resulting electric field distribution in the vicinity of a layered elliptical 
cylinder with inner (outer) permittivity e\ — 2(e2 = 3) and (b) with inner (outer) permittivity 
ei = 3(e2 = 2). The light is incident from the — x direction, in the x — y plane, normally to the z 
axis. 
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